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We perform a detailed numerical investigation of the dynamics of a single component broken 
symmetry A<^'* field theory in 1+1 dimensions using a Schwinger-Dyson equation truncation scheme 
based on ignoring vertex corrections. In an earlier paper, we called this the bare vertex approxima- 
tion (BVA). We assume the initial state is described by a Gaussian density matrix peaked around 
some non-zero value of (0(0)), and characterized by a single particle Bose-Einstein distribution 
function at a given temperature. We compute the evolution of the system using three different 
approximations: Hartree, BVA and a related 2PI-1/N expansion, as a function of coupling strength 
and initial temperature. In the Hartree approximation, the static phase diagram shows that there 
is a first order phase transition for this system. As we change the initial starting temperature of 
the system, we find that the BVA relaxes to a new final temperature and exhibits a second order 
phase transition. We find that the average fields thermalize for arbitrary initial conditions in the 
BVA, unlike the behavior exhibited by the Hartree approximation, and we illustrate how {4>{t)) and 
depend on the initial temperature and on the coupling constant. We find that the 2PI-1/N 
expansion gives dramatically different results for {(l>{t)). 

PACS numbers: 11.15.Pg,03.65.-w,11.30.Qc,25.75.-q 



I. INTRODUCTION 

Recently there has been much effort in finding approx- 
imation schemes to study the dynamics of phase tran- 
sitions that go beyond a leading order in large-N mean 
field theory approach. This is an important endeavor 
if one wants a first principles understanding of the dy- 
namics of quantum phase transitions. In a previous set 
of papers, we studied in quantum mechanics jl] , as well 
as 1-1-1 dimensional field theory |^, the validity of a 1/iV 
motivated resummation scheme, which we called the bare 
vertex approximation (BVA). 

The long-term goal of this work is finding approxima- 
tion schemes which are accurate at the small values of 
N relevant to the case of realistic quantum field theo- 
ries: N=4 for the linear sigma model, and N=2 for the 
Walecka model. In order to maximize the possible differ- 
ences between approximation schemes we choose to study 
the 0(N) model for iV = 1. Based on our previous studies 
of the quantum mechanical version of this model Q we 
expect that by increasing N the differences will dimin- 
ish. This is due to the fact that the Schwinger-Dyson 
formalism is related, but not identical, with approxima- 
tions based on the large-N expansion. 

This paper presents the first quantum-mechanical dy- 
namical calculation in 1-1-1 dimensions, using the BVA, 
of the explicitly broken symmetry case, where the order 
parameter at i = is non-zero. In 1+1 dimensions, it 
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is known that there is no phase transition in this model 
except at zero temperature On the other hand, in 
two dimensional systems having Berezinski-Kosterlitz- 
Thouless type transitions ^, ^ , the large-N expansion 
can give qualitatively good understanding of the correla- 
tion functions, even when it gives the wrong phase tran- 
sition behavior |^]. As the dimensions increase, the mean 
field critical behavior becomes exact in four dimensions 
and thus we expect that the approximation presented 
here should improve as we increase the dimensionality. 
Thus we should think of the model used here as a "toy" 
model for demonstrating some of the features expected 
to be true in 3+1 dimensions such as the restoration of 
symmetry breakdown at high temperatures and thermal- 
ization of correlation functions. The model we are ul- 
timately interested in is the linear sigma model in 3+1 
dimensions for the case of broken symmetry at finite tem- 
perature. This model we studied earlier using a large-N 
approximation ||, Since the BVA contains scat- 
tering contributions, it cures many of the problems as- 
sociated with mean field methods. Here we are able to 
follow the evolution of the system through a phase tran- 
sition and study the thermalization of the system. We 
illustrate in this paper the results of such calculations. 

A parallel set of investigations by Berges et. al. |l^, 
[l^ , |l3| , |l^ , |l5|, have looked at a related approximation 
based on the two-particle irreducible expansion (which 
they call 2PI-1/N). These investigators have pointed out 
that when there is broken symmetry, the BVA contains 
terms not included in the 1/iV resummation at next to 
leading order. In this paper, we present the first quantum 
calculations which compare the BVA with the 2PI-1/N 
for the broken symmetry case. Recently we were able to 
show that for a classical finite temperature X(j)'^ field the- 
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ory in 1+1 dimensions, the BVA gave a better description 
of the time evolution of {(p{t)) than the 2PI-1/N expan- 
sion and provided excellent agreement with exact Monte 
Carlo simulations H, |l^ . 

In this paper we look at quantum evolutions in 1+1 di- 
mensions, starting with a Gaussian density matrix, and 
study how the evolution of {(j>{t)) depends on the initial 
conditions and the value of the coupling constant. In the 
classical domain, the coupling constant dependence can 
be scaled out, which is not possible in the quantum case 
we consider here. Since we have not determined the effec- 
tive potential in the BVA approximation, we rely on the 
Hartree approximation effective potential to guide our 
study. The Hartree potential however, indicates that the 
system should undergo a first order phase transition. In 
addition, in the Hartree approximation, the fields never 
thermalize. We find that the BVA cures these serious 
problems. In this paper, we show evolutions of the sys- 
tem as a function of the initial "temperature" parameter 
and the coupling constant, and, since in the BVA, the 
fields thermalize, we can follow the system through what 
appears to be a second order phase transition. We also 
compare our results with those of the 2PI-1/N expansion 
in the quantum domain, and find that the two approxi- 
mations diverge dramatically for the behavior of {(j){t)). 

We derive the BVA equations for the general iV- 
component A [ 0f (a;) ]^ field theory in Sections |l| and [H. 
We then specialize to the case = 1, and in Section IV, 
we derive the Hartree approximation phase diagram. In 
Section 0, we discuss the initial conditions we choose for 
this problem. Numerical results are shown in Section [v|, 
and conclusions discussed in Section VII. 



and the constraint ("gap") equation for x{^) 



A 



X{x) = + 2]v" 'l^i{x)4>i(x) . 



(4) 



Throughout this paper, we use the Einstein summation 
convention for repeated indices. 

The BVA truncation scheme of the Schwinger-Dyson 
equations is most easily obtained from the 2PI effective 
action |l^, |l9|. Other approaches leading to these 
equations are found in ||l|, |l^ . Using the extended fields 
notation, = [x{x), (t^iix), (j)2{x), . . . , 4>n{x)] , the ef- 

fective action for the evolution can be written as: 



1 + -Trln[G-i] 



(5) 



where T2[G] is the generating functional of the 2-PI 
graphs, and the classical action in Minkowski space is 



(fx 



]^(t),{x)[U + x{x)]^^{x) 



2ff 



— X{x) 
9 



(6) 



Here and in what follows we let g = X/N. 

The integrals and delta functions Sc{x,x') are defined 
on the closed time path (CTP) contour, which incorpo- 
rates the initial value boundary condition ||2^, ^ ^ . 
The approximations we are studying include only the 
two- loop contributions to 

The Green function G^^^[0](a;, x') is defined as follows: 



II. THE CLASSICAL ACTION AND TIME 
EVOLUTION IN THE BVA 

The classical action for X(j>'^ with N fields (z = 1 . . . N) 



IS 



1 r 



8N 



2A 



(1) 



For the purposes of our resummation scheme which is 
motivated by 1/A'^ considerations it is useful to consider 
the alternative action 



2 L-'^ 



-^ + ^^ x{x) 



(2) 



which leads to the Heisenberg equations of motion 

[D+x{x)](b,{x) = 0, (3) 



Oa/3l 



[X, X 



5^Sc 



S(t)aix) Scjjpix') 



D^\x,x') K~Hx,x') 
K-l{x,x') G-l{x,x') 



(7) 



where 



-Do ^{x,x') = -g5c{x,x') , 
Goy- [x] {x, x')^[U+ x{x) ] 5^.J5c{x, x') , 
K-lmx,x') = k-l[4>]{x,x') = c^,{x)5c{x,x') . 

The exact Green function Gaf3[j]{x,x') is defined by: 

' D{x,x') Kj{x,x'y 



G^f3[j]ix,x') - (^^(^^^z) gUx,x') 



The exact equations following from the effective action 
Eq. (^), are: 

[D + x{x)]Mx)+K,{x,x)/i = 0, (8) 
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and 



where 



G^^{x, x') — Gq^p{x, x') + T,ai3{x, x') 



(9) 



J:a0{x,x') 



2 ST2[G] 
i 6Gap{x,x') 



(10) 



In the BVA, we keep in r2[G] only the graphs shown in 
Fig. which is exphcitly 



dxdy 



Gtj{x,y)Gj,{y,x)D{x,y) (11) 
+ 2Ki{x,y)Kj{x,y)G.ij{x,y) . 
The self-energy, given in Eq. (px|), then reduces to: 

Il{x, x') = — Gmn{x, x') Gmn{x, x') , (12) 

(x, X ) — % Km {x^ X ) Girii (x, X ) , 
flj(x, X ) — t Km{x^ X ) Grniix^ X ) , 

(x, x') — i[ Gij (x, x') D{x, x') + Ki(x, x') Kj{x, x') ] . 

As discussed in detail in Ref. the second graph 
in Fig. is proportional to 1/A^^ and is ignored in the 
2PI-1/N expansion. Our recent simulations in the clas- 
sical domain showed that the BVA gave a more accurate 
determination of {(j){t)), and we will concentrate in this 



paper on the BVA except to point out with explicit re- 
sults that in the quantum domain the differences between 
the BVA and the 2PI-1/N expansion grow with increas- 
ing coupling constant g (for the case N — 1 studied here) 
and that unlike the BVA, the 2PI-1/N expansion does not 
track the average of the Hartree result. 



+ 
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FIG. 1: Graphs included in the 2PI effective action 



III. UPDATE EQUATIONS FOR THE GREEN 
FUNCTIONS 



We notice from the definitions of the matrices repre- 
senting Gai3{x,x') and G~p{x,x'), that the matrix ele- 
ments are not inverses of one another, but instead satisfy 
schematically: 



D-^D + K^^Kk = Sc, 
K'D + GT^'K, = 
D-'K,+K^'Gk,=0, 
^i^Kj + G-^Gkj = S^jSc . 
Inverting Eq. (^, we find: 



(13) 



with 



D(x,x') = —g5c{x,x)+g I dxi n'(x, xi) D(xi, x') 



Gy(x,x') = Goij(x,x')(5y - / dxi / dx2 Goifc(a;, xi) S'^; (xi, X2) (x2, x') , 

JC Jc 



Ki{x,x) 



dxi / dx2 
c Jc 



(x,xi) K^l + Vlk (xi, X2) Gfci(x2, x') 



n'(x,x') = n(x,x')— / dxi / dx2 
Jc Jc 

E-fe(x,x') = i;ife(x,x')- / dxi / dx2 

Jc Jc 



(x,xi) 



^0 ki + 



{X1,X2) 



(x,Xi) Do^+n (xi,X2) 

r 



(X2,x') , 
(X2,x') 



(14) 
(15) 
(16) 

(17) 
(18) 



These update equations must be solved in conjunction 
with the one-point functions, Eqs. (H). 

For a practical implementation of the above approach 
we need to solve for D2{x,x') and G2.ij(x,x'), the in- 



verses of 



^0 ' + n 



(x,x') and [G^ ^ + S^j (x,x'), re- 



spectively. We have 



D2(x,x') = - g5c{x,x') 



+ 9 J dxin(x, xi)I?2(a;i, a;') , 



(19) 



G2,ij{x,x') = Goij{x,x') (20) Thus we obtain the equations of motion 

- dXi dX2Goik{x,Xi)i:kl{xi,X2)G2,lj{X2,x') . 

Jc Jc ( 1 _ 

I [ □ + x{x) ]5ik + g Gki{x, x)/ij(t)k{x) + Ki{x, x)/i = 0, 



We also perform the following substitutions 

D{x,x') = -gSc{x,x') + D{x,x'), (21) 
D2{x,x') = -g6cix,x') + D2{x,x') , 
Ki{x,x') = g4>k{x)Gki{x,x') + Ki{x,x') . 



X{x) = -M^ + ^<1>K^) + I [Guix,x)/i] . (22) 



and the update equations for the Green functions 



and 



with 



D(x,x') = —g^ll'{x,xi)+g J dxi Il'{x,xi)D{xi,x') , 
D2{x,x') = -g'^U{x,xi) + g j dxi'n.{x,xi)D2{xi,x') , 
Gij{x,x') = Goij{x,x') - g j dxi Goikix,xi) (j)k{xi)(j)i{xi) + Gki{xi,xi)/i Gij{xi,x') 

- dxi dX2 Goik{x,Xi)f:'ki{xi,X2)Glj{X2,x') , 

Jc Jc 

G2ij{x,x') = GQij{x,x') - g j dxi Goik{x,xi)[G2ki{xi,xi)/i]G2ij{xi,x') 

dxi / dX2 Goik{x,Xi)f:kl{xi,X2)Glj{X2,x') , 



Ki{x,x') = - j dxi D2{x,xi)<j)k{xi)Gki{xi,x') 

+g / dxi Qkix,xi)Gki{xi,x') - / dxi / dx2 D2{x,xi)^lk{xi,X2)Gki{x2,x') , 
Jc Jc Jc 

U'{x,x') = IL{x,x') - 4)k{x)G2,kl{x,x')(j)i{x') - I dXi j dX2 ^k{x,Xi)G2,kl{Xl,X2)^l{X2,x') 

Jc Jc 

- I dxi (l)k{x)G2,ki{x,xi)Qi{xi,x') - / dxi Q.k{x,xi)G2,ki{xi,x')(j)i{x') , 
Jc Jc 



IC JC 

Ejfe(x', x') = i [Gik{x, x) D{x, x') + Ki{x, x) Kk{x, x')] 



^k{x,x') = Ysik{x,x')- 4>^{x)D2{x,x')(t)k{x 



')- I da^i / 
Jc Jc 



dX2 Q.i{x,X-i)D2{x-i_,X2)^k{X2,x') 



■ I dxi 4>i{x)D2{x,xi)Q.k{xi,x') - I dxi Q.i{x,xi)D2{xi,x')(j)k{x') 
Jc Jc 



+g 



4>i{x)ilkix,x') + ni{x,x')4>k{x') +g I dxi xi) rife(xi, x') . 

c 



(23) 
(24) 
(25) 

(26) 



(27) 



(28) 

(29) 
(30) 



For computational purposes, it is suitable to make one more transformation of the 0,, and G2ij equations. We 
write the equivalent integro- differential equation for Gij as: 

|[a + x{x)]Sik + g[(t>i{x)(t>k{x) + Gik{x,x)/i]^Gkj{x,x') = Sc{x,x')Sij - j dxif:ik{x,X2) Gkj{x2,x') . (31) 



We specialize now to the case N = 1. It is convenient 
then to introduce the following equations 

[D+xiix)]cP{x)+K{x,x)/i = 0, (32) 



X,{x)=-i,^ + ^cl^\x)+'^Gix,x)/i, (33) 
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together with redefinitions for Go{x, x') to work with the 
equations for G{x,x') and G2{x,x'), respectively. We 
have 

[a + x2{x)]Goix,x')^6c{x,x') (34) 
[a + xiix)]Go{x,x') =Scix,x') , (35) 



with 



X2{x) = -^^^ + ^U'{x) + G{x,x)/^] , (36) 



Finally, the modified equations are given by: 

G{x,x') = Go{x,x') (37) 



dxi / dx2Go{x,xi)Y:'{xi,X2)G{x2,x') , 
c Jc 



G2{x, x') — Go{x, x') 



Axi I Ax2Gq{x,Xi)T.{xi,X2)G2{x2,x') 

c Jc 



(38) 



IV. HARTREE PHASE DIAGRAM 

It would be useful to have available the effective po- 
tential for the BVA approximation from Eq. (|^) to use 
as a guide for starting out the BVA solutions. However 
solving the self-consistent equations of the BVA and con- 
structing the thermal effective potential is a formidable 
task, and has only been recently considered for the 
simpler loop approximation to for the single field 
{N = 1) case (see Refs. |5|). Therefore in this pa- 
per, we consider the effective potential for the simpler 
Hartree approximation for a single field. 

The effective action in the Hartree approximation can 
be written in the form 



d2a;|-l0[n + x]0+^0' (39) 



This action gives the Hartree equations of motion: 

[a + x{x)-X.p^{x)](l){x)^0, (40) 
X(:^:) = V + Y'/''(2;) + YTr[GoA], 

with Gq^{x,x') = [n + xix)]d{x,x'). 

The effective potential for this action is given by: 

/•+00 dfc r 2 1 
+ 2^ r ^' + /? ^''[^ " cxp(-/3c^fe)] I , 



where cuk 
dx 



^Jk?' + X- We note that the requirement: 
+°° dfc 2nu + 1 



Y-3a(^ + '^) + 



27r 2a;fc 



= 0. 



leads to the gap equation: 
9 3A 2 



X 



2^ 



1 



/o ^" ^fe 

where 71^ = 1/ [e^'^'' — 1] . The above equations are infinite, 
so to renormalize them, we introduce a cutoff at fc = ±A, 
and introduce a quantity > 0, defined by: 



/i 



3A 

T 

3A 

in 



dk 



1 



27r Vfc2" 
ln( A/m) . 



(41) 



Recall that /i^ > 0. Subtracting (|4^) from the gap equa- 
tion gives: 

9 3^ ,2 

'"^ dfc / 2nfc + 1 1 



X = - m 



(42) 



3A 
2 



2^ V yfc^T^ 



3A 
T 



A 



— In(mVx) 



dfc rife 

TT Wfe 



which is now finite. The Hartree potential, is renormal- 
ized at T = by first renormalizing the partial derivative: 



dVii[(b,x] 
dx 



2 



1 

3A 



(x- 



dfc 
2^ 



2v^ 



(y - 

2 3A 



X 2Vfc2 + 

-!- In(mVx) , 
an 



where we have used ( |4l| ) to make the equation finite. 
Partially integrating, we obtain the renormalized Hartree 
effective potential: 

^h[0,x] = ix'/''-^0'-^(Y+m2x) (43) 
1 

+ 8^^ 



X~X ln(m^/x) 
1"°° dk 1 

/ ^ ln[l-exp(-/3c.fe)], 



where we have added back in the finite temperature- 
dependent part. This equation is to be solved with x 
satisfying the renormalized gap equation (p^). In prac- 
tice, it is useful to solve both of these equations paramet- 
rically as a function of x- 

The physical (renormalized) mass is given by the sec- 
ond derivative of the effective potential, evaluated at the 
minimum. The minimum occurs at 



= v{x-W) = 0. 
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This implies that for the symmetry-breaking solution, 
X = Aw^. Thus, from the gap equation (^, the po- 
sition of the minimum and the mass parameter are 
related by: 



A 



3A 
T 



dfc nfe(Au^)"' 



(44) 



The renormalized mass is defined by: 



X - 3A02 + 











d(/)2 








4)—v 



^From the gap equation (|4^), we find: 
dx 1 3Aw 



1 + 3[1 + /(A«2)]/(8^Z;2) 

where /(Aw^) is the finite integral: 



/ X = X / dfc 3 



Thus the renormalized mass can be computed from: 
^;2-3[l + /(A«2)]/(4^) 



m 



R 



A 



l + 3[l-H/(Aw2)]/(87rw2 



The critical temperature Tcr is defined by the simultane- 
ous solutions of: 



An 



and Eq. 



At T = 0, we notice that unless v > 



3/(47r), one cannot have a symmetry- breaking solution 
in this approximation. The effective potential as a func- 
tion of temperature T can be computed numerically. In 
Fig. H we depict the effective potential dependence on 
the coupling constant at a fixed temperature, T = 0.1. 
In Fig. ^ we fix the coupling constant at A = 7.3, and 
depict the dependence of the effective potential on the 
temperature. This particular value of A was used in our 
study of the dynamics of disoriented chiral condensates 
(DCC) in 3 + 1 dimensions in the leading order in large-N 
approximation ^, [l^. Here, the phase transition occurs 
with Tcr ~ 0.878. We see that the phase transition is 
first order with the vacuum value v w 0.635. The value 
of X Ett this vacuum is Aw^ « 2.94. 



V. INITIAL CONDITIONS 

For the purpose of this study, we will assume that ini- 
tially (at t = 0) the system is described by a Gaussian 
density matrix. Thus initially the field equation and two- 
point equation is that of the Hartree approximation. The 
field equation obeys (pO|), with Go{x,x') satisfying: 



(45) 




0.2 0.4 0.6 0. 



1 1.2 1.4 1.6 



FIG. 2: We set T = 1, and plot the Hartree effective potential 
as a function of the coupHng constant. 




FIG. 3: We set A = 7.3, and plot the Hartree effective poten- 
tial as a function of temperature. 



We can solve this Green function equation by introduc- 
ing a set of quantum fields (f>o{x), satisfying canonical 
commutation relations [(j)Q{x), (fioix') ] = i S{x — x') , and 
obeying the homogeneous differential equation 



[a + x{t)]M^)^0- 

In terms of these fields we have 



(46) 



Go{x,x') = i{Tc{M^)M^')}) , 

= Gyix,x')ecit,t') + G<ix,x')ec{t',t) , 

+ 00 Ah, 

— Goik;t,t')e'''^''~''K (47) 
-oo 27r 

We next expand these operators in Fourier mode func- 
tions 



dfc 



(48) 

where the mode functions /j, (t) satisfy: 
[df+u;l{t)]Mt)^0, LOu{t)^^B+x{t), (49) 
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and the Wronskian condition: f^{t) f^. (<) — f^{t) (t) = 
—i. The operators and ajj, satisfy the usual commu- 
tation relations: [oq J,, ajj,, ] = 27r(5(fc — fc'). We will take 
our initial density matrix such that: 

(flofcaofc' ) = nk2Ti 5{k - k') , {a^^a^^,) = 

i%kalk') = K + l)2TT5{k-k') , (aj^ajfc, ) =0, 

(50) 

where nfc = l/[exp(/3oWfc(0)) — 1 ]. Here Tq = l/Po is just 
a parameter for the initial Gaussian density distribution, 
and is not the true temperature of the interacting system. 
In fact, the system will not be in equilibrium at t = 0, 
but will arrive at a final temperature T after it has come 
to equilibrium. 

Solutions for the mode functions /j, (t) arc of the form: 



fkit) 



-i J*nkit')dt' 



(51) 



where Vlk (t) satisfies the non-linear differential equation 



+ nUt)=ujlit). (52) 



The first order WKB solution for {t) is then given by 
^lk{t) = 0Jk{t). We take these solutions for our initial 
conditions, so that at i = 0, 




f]fc(o) = c^fc(o) = Vfc2 + x(o), 

f2fc(0)=^^fe(0)=x(0)/2a;fc(0). 
This means that 



/fe(0) = 1/V2c^fe(0) 
^fe(O) 



/fc(0) 



2cjfc(0) 



icjfe(O) 



(53) 



(54) 



Wc still need to find the value of Xo(0)- This is given by 
the Hartree self-consistent solutions of: 



X(0) 



3A 
3A 



^(0) 



3A 
T 



dfc 2nk 4- 1 



^0) 



2^ V^^+X(0) ' 
(55) 



3A 
Y 



dfc 
2^ 



2nfc 



Where we have used Eq. (|4l|). 

So, for our case, Fourier transforms of the Green func- 
tions at i = are given by: 

Go>(fc;t,OA - f kit) flit') {nu + 1) + /,*(t)/fe (f) "fe , 

Go<(fc;t,OA = /fc(t)A.(0 (rife + 1) + A.(0/fc(i')"fc • 

(56) 

These results, together with Eq. (|5^), determine the val- 
ues of Go >(fc; i, t'), and all it's derivatives, &tt = t' = Q. 



VI. NUMERICAL RESULTS 

We choose initial conditions for the two-point functions 
as described in the last section. In all our simulations we 
set the renormalized mass parameter m, as defined in 
Eq. (^Tj) , to unity. We have verified numerically that our 
results are independent of the cut-off parameter A for 
values of A between Stt and 47r. The numerical proce- 
dure for solving the BVA equations is described in detail 
in Rcfs. ||2^, |2^. The energy is numerically conserved to 
better than five significant figures for all our simulations. 
The calculations are carried out entirely in momentum 
space, and the results are free of artifacts related to the 
finite volume of a lattice in coordinate space. For a dis- 
cussion of differences between the continuum and the pe- 
riodic lattice approach, respectively, we refer the reader 
to our previous paper, Ref. pst . 

We start by choosing an initial temperature Tq — 0.1 
so as to bring out the quantum nature of the dynamics. 
The Hartree effective potential for this value of Tq is a 
slowly varying function of the coupling constant, with the 
minimum value of at large A being between 0.6 and 1.0 
(see Figs. || and Thus we want to choose small initial 
values of 0(0) = 0.4 and take 7r(0) = 4>{Q) = so that 
we start below the height of the barrier. Then we should 
see (f) moving to the opposite side of the well and then 
settling down at the potential minimum position. We 
show the results of this calculation in Figs. ^ and || for 
several values of the coupling constant. We notice that 
the position of the BVA minimum is located between 0.55 
and 0.60 instead of the Hartree average value of 0.8. We 
also see that the final value of x does seem to be linear 
in A as in the Hartree result. However, the ratio xli^'^'^) 
is one for the Hartree approximation, but is about 2.4, 
with a less than 1% error, for the BVA simulations. 




FIG. 4: Coupling constant dependence: We set To = 0.1, and 
plot {(j>{t)) for various values of A. 

Next we set A = 7.3, which is the phenomenological 
choice for the the linear sigma model, with (/)(0) — 0.4 
and 7r(0) = 0(0) — 0, and study the dependence of {(j){t)) 
and {x{t)) as a function of the initial temperature Tq. 
The results are shown in Figs. ^ and ^. We note im- 
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FIG. 5: Coupling constant dependence: We set To = 0.1, and 
plot (xii)) for various values of A. 
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FIG. 7: Temperature dependence: We set A — 7.3, and plot 
(x(i)) for various values of Tq. 



mediately from Fig. ^ that for small values of Tq, {(t>{t)) 
equilibrates to non-zero values, but for large values of Tq, 
{4>{t)) equilibrates to zero, as expected from the Hartree 
effective potential. Fig. |^ shows that x(t) equilibrates to 
different values which depend on Tq. 
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FIG. 8: We are increasing the kinetic energy: We set A — 7.3, 
To = 0.1, and plot {(j>{t)) . 



shown in Figs. ^ and 



FIG. 6: Temperature dependence: We set A — 7.3, and plot 
{4>{t)) for various values of Tq. 

Since we have chosen the case 7r(0) = 0(0) = 0, 
with (/)(0) just under the barrier height, the transition 
to the equilibration point is very slow, i.e. equilibration 
is reached without too many exciting features. It is in- 
teresting to give the system a little initial kinetic energy 
so that structure is introduced in the dynamics but the 
equilibration value of {(j){t)) remains the same. We illus- 
trate this in Fig. |[ 

The plots of the BVA order parameter in Fig. ^ indicate 
that for very low values of Tq, the order parameter {4>{t)) 
approaches a non-zero constant. For very large values of 
To, the order parameter goes to zero, as expected. Some- 
where between Tq — 1.0 and Tq — 2.5, there seems to be 
a phase transition. In order to study this in more detail, 
we have carried out BVA simulations for temperatures 
between Tq = 1.5 and Tq = 2.5 at 0.1 intervals. The re- 
sults for the order parameter and its first derivative are 




To=l-5 
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To=1.8 



To=2.1 
To=2.2 
To=2.3 
To=2.4 
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FIG. 9: Detailed temperature dependence: We set A = 7.3, 
and plot {4>it)) ■ 

The equilibration points for the order parameter in 
Fig. ^ are at almost equally spaced intervals for Tq < 2.4, 
and there is no evidence of a jump in the value of the or- 
der parameter from a finite value to zero, as would be 
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FIG. 10: Detailed temperature dependence: We set A — 7.3, 
and plot {(j){t)) . 



FIG. 12: Comparison of various methods: We set A — 1.0 and 
To = 0.1, and plot (xW)- 



required by a first order phase transition 
((^(t)), shown in Fig 
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The plot of 
reveals the complex nature of 
the dynamics: {4>{t)) oscillates and gradually approaches 
zero, when equilibration occurs. Oscillations around a 
zero value of suggest the trapping in a potential 

well, while the wiggles at negative values of may 
be due to changes in the position of the local minimum of 
the effective potential as a function of temperature and 
time. Since the chosen initial conditions do not represent 
an equilibrium state for the interacting system, the dy- 
namics toward a final state of equilibrium is accompanied 
by changes in the effective temperature. In other words, 
the effective potential is not such a good guide to the dy- 
namics. In analyzing these figures, it becomes apparent 
that a) there is no abrupt change that would support the 
first-order labeling of the phase transition, and b) there is 
a qualitative change happening near Tq = 2.4, where the 
plot begins crossing the lower temperature curves, 
and the local oscillations in completely disappear. 

Therefore, we conclude that the phase transition, as cal- 
culated using the BVA, is probably not first order. 
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FIG. 13: Comparison of various methods: We set A : 
To = 0.1, and plot {<?!>(*)). 
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FIG. 14: Comparison of various methods: We set A = 7.3 and 
To = 0.1, and plot {x(i))- 
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FIG. 11: Comparison of various methods: We set A = 1.0 and 
To = 0.1, and plot {(l>{t)). 

It is known that this model in two-dimensions has 



no phase transition at finite temperature p9| . At 
zero temperature, there is a second order phase tran- 
sition p^pll Explicit Monte Carlo lattice calcu- 
lations T33,^3|have shown that indeed theory in 
two-dimensions and at zero temperature is non-trivial 
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at least when the continuum hmit is reached from the 
broken symmetry phase, and that the symmetry is fully 
restored at high temperature. It is known however that 
approximate lattice calculations (such as the variational- 
cumulant expansion method [ ^5| , p^), which are designed 
to study scalar 0^ theory in 3+1 dimensions on the lat- 
tice, may erroneously indicate the presence of a second- 
order phase transition at finite temperature in 1-1-1 di- 
mensions. This is probably due to the fact that the ex- 
pansion is only carried out to third order. Much in the 
same way, the Hartree approximation exhibits a first or- 
der phase transition with T^r « 0.878. When going be- 
yond the Hartree approximation (using the BVA) we still 
find a phase transition, but the BVA relaxes the order 
of the phase transition. At very low temperature, we do 
(apparently) find a non-zero value of the order parameter 
as t — > oo, which is the exact result at zero temperature. 
So we conclude that the BVA at low temperature is see- 
ing effects that might occur in higher dimensions, even 
though it is not technically correct in 1+1 dimensions. 

We next compare three different approximation meth- 
ods: Hartree, the 2PI-1/N expansion, and BVA. Here we 
choose a very low initial temperature of Tg = 0.1 in order 
to emphasize quantum effects in the dynamics. Again, we 
start with ((/)(0)) = 0.4 and MO)) = 0. We show results 
for A = 1 in Figs. |ll| and y and A = 7.3 in Figs. |l| 
and |l^ As expected, we find that the Hartree approxi- 
mation leads to oscillation about the Hartree minimum 
without equilibration. The BVA results track the Hartree 
curve, except with damping, and go to a non-zero value 
as t — + cxD, which is the exact result at zero temperature. 
The 2PI-1/N expansion goes to a zero value of the or- 
der parameter, in agreement with the exact result of no 
phase transition for finite temperature. Results for {x{t)) 
are similar for both approximations, in agreement with 
expectations based on our experience with the classical 
limit of these approximations |Lq | . The disagreement be- 
tween 2PI-1/N and BVA is more pronounced at larger 
values of A, as shown in Figs. n3 and 04. 



where we set 



A{t) 



(59) 



with m^jf(i) the effective mass and Teff(t), the effective 
temperature at time t. The factor A{t) comes from a 
wave function renormalization at each time t. So we can 
determine nk(t) and ujk(t) from the relations: 



UJk{t) 

nk{t) = 



d^G{t, t'-k) 



dt dt' 



d^G{t, t'-k) 



dtdt' 



G{t,t;k) 



G{t,t;k) 



(60) 
(61) 



^From (pOj), we notice that uJk is a ratio of Green func- 
tions, and thus any (finite) wave function renormaliza- 
tion will cancel. However from (^Tj), Uk is directly pro- 
portional to the Green function, which will have a finite 
wave function renormalization when restricted to the sin- 
gle particle contribution. 
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FIG. 15: For A 

of LOk- 



7.3 and To = 2.5 we plot the time evolution 



VII. THERMALIZATION 

The BVA leads to thermalization of the system. In 
order to have a measure of the semi-static thermody- 
namic properties of the system, it is reasonable to fit 
our time-dependent Green functions to those appropri- 
ate to a free field with frequency ojk (t) and Bose-Einstein 
distribution function nk{t), see also That is, we 

equate the Fourier transform of the BVA Green func- 
tions Gk{t,t) and d^Gk{t,t')/dtdt' to the corresponding 
free-field cases: 



Gk{t,t) = [2nfe(t) + l] 



1 



d^Gk{t,t') 



dtdt' 



[2nk{t) + 1] 



2ujk{t) ' 

UJk{t) 



(57) 
(58) 



First let us concentrate on LL>k{t). We show in Fig. 
a plot of uik{t) as a function of k for various values of t, 
as calculated from Eq. (|60|), for A = 7.3 and Tq = 2.5, 
which is close to the critical temperature. We fit the 
LOk{t) data to determine the values of rn^Q^t). These re- 
sults are compared to the value of x(i) in Fig. |l^. We 
see here that self-energy corrections to the effective mass 
reduce the effective mass to about 25-40% of x{t)- Fur- 
ther calculations for different initial temperatures allow 
us to conclude that the correction slowly decreases with 
increasing initial temperature. 

Now that we have the effective uJk we can find out if 
the particle number density nk(i) has the simple Bose- 
Einstein form given in Eq. (p9[). We use a nonlinear 
fitting procedure to obtain the parameters T'eff(t) and 
A{t) from the data generated by Eq. (|6l|). The results 
are shown in Fig. |l^ for A = 7.3 and Tq = 2.5. In 
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FIG. 16: For A = 7.3 and To = 2.5 we compare m'^{t) and 

xit). 



FIG. 18: For A = 7.3 and To = 2.5 we plot the time evolution 
of rifc, rcn, after renormalization. 
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FIG. 17: For A = 7.3 and To = 2.5 we plot the time evolution 
of rifc, before renormalization. 



FIG. 19: For A = 7.3 and To = 2.5 we plot the time evolution 
of the temperature Tk,Tan{t)- 



Fig. |l9| we show the renormalized densities, given by 
nYcnk{t) — nk{t)/A{t). Notice that nk(t) is very simi- 
lar to a Bose-Einstein distribution which starts out at a 
temperature of TcS = To, increases in amplitude due to a 
temperature increase, then falls back to a lower temper- 
ature, but larger mass. 

In Fig. |9| we show the effective temperature profile, 
Tfc,rcn(t) inferred from nk,venit) 

Trcnk{t) = LOk{t) / In [ ( 

^rcn k W + l)/KcnfcW)] . (62) 

We notice that at short times, T^cnkit) has to readjust 
to the effects of the non-Gaussian corrections, but then 
it settles down to a new temperature which is relatively 
independent of the momentum k. We also notice that 
we do have quite a bit of particle production: after an 
initial spike, the particle density number relaxes to an 
equilibrium value. 

In Fig. ^ we show the dependence of Teff on the initial 
temperature Tg, for fixed A = 7.3, and in Fig. the 
dependence of Teg on the coupling constant A for fixed 
To = 0.1. 




\— — , , , , , , , , 1 
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FIG. 20: For A = 7.3 we plot the time evolution of T^«(t) for 
various initial temperatures Tq. 



VIII. CONCLUSIONS 

In conclusion, we have shown that the BVA indeed 
leads to equilibration of the one- and two-point functions 
and thus remedies this deficiency of the Hartree approxi- 
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FIG. 21: For To = 0.1 we plot the time evolution of Tcff (t) for 
various values of the coupling constant A. 

mation. The nature of the phase transition changes from 
first to second order when going from Hartree to BVA 
which is the correct behavior (as a function of couphng 
constant) at zero temperature but is incorrect at higher 
temperature. However, such a phase transition at finite 
temperature is expected in the 3+1 dimensional case that 
we are interested in modeling in the future. The 2PI-1/N 
approximation does not see this phase transition (which 
is not present in the exact theory) and does not track the 
average of the Hartree result. In higher dimensions, one 
expects the Hartree result for the order parameter to be 
correct on the average but not have the property of equi- 
librating. So deciding which approximation is more phys- 
ical will have to wait for a 3-1-1 dimensional simulation 
to see if both approximations show a phase transition. 

In the classical regime, where we could do an ex- 
act Monte Carlo calculation, we found that the BVA 
works better than the 2PI-1/N expansion in capturing 
the dynamics of {(j>{t)), but that the difference was not 
very great. In the quantum domain, however, these 
two approximations diverge at low temperature with the 
strength of the coupling constant, and the 2PI-1/N ap- 
proximation no longer tracks the average of the results of 
the Hartree approximation. Since we do not have exact 
calculations in the quantum regime, we cannot make any 
strong conclusions about this divergence. 

It is interesting to note that technical issues related 
to the nature of the phase transitions in 1-1-1 dimensions 
are not uncommon. There is a previous example of just 
this very type of reasoning being successful. In trying 
to understand the QCD chiral phase transition at finite 



chemical potential and temperature, it was important 
to have a model where the phase transition mimicked 
what is known in 3+1 dimensions. Such a model was 
a four-fermi model (the Gross-Neveu model) in 1+1 di- 
mensions in leading order in large-N, which had a similar 
phase structure to two-flavor QCD. This model was then 
used as a testing ground for studying the effect of first 
and second order phase transitions, with a critical point, 
and showed some qualitative differences between the two 
types of transitions (see Ref. 0). Once again, the exact 
1+1 dimensional model does not have a phase transition. 
But this warmup problem then allowed us to go to 2+1 
dimensions p8| , where it is known that the leading order 
large-N approximation and the exact theory have simi- 
lar phase diagrams, as verified by lattice simulations in 
2+1 dimensions We submit that this is the way 

these 1+1 dimensional results should be understood: as 
a model having certain properties when treated in this 
approximation and as a testing ground for codes which 
will then be generalized to higher dimensions where it is 
expected the approximation will correspond to the known 
behavior. The main purpose of the quantum simulations 
we present here is to get experience in getting codes work- 
ing in lower dimensions that qualitatively do what we 
expect to see in 3+1 dimensions. 

As a result of the simulations presented here, we are 
confident that we can now study the chiral phase tran- 
sition in 3+1 dimensions in the linear sigma model and 
describe the competition between the expansion of the 
plasma and the equilibration tendencies. This will al- 
low us to see whether some of the phenomena present 
in Hartree (and/or large-N) approximation, such as pro- 
duction of disoriented chiral condensates and distortion 
of pion and dilepton spectra, are still present in spite of 
the forces that lead to thermalization. 
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